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ABSTRACT 

We study orbit of a pressure-confined cloud in the broad-line region (BLR) of active 
galactic nuclei (AGNs) when the combined effects of the central gravity and anisotropic 
radiation pressure and the drag force are considered. Physical properties of the inter¬ 
cloud gas such as its pressure and dynamic viscosity are defined as power-law functions 
of the radial distance. For a drag force proportional to the relative velocity of a cloud 
and the background gas, a detailed analysis of the orbits is performed for different 
values of the input parameters. We also present analytical solutions for a situation 
where the intercloud pressure is uniform and the viscosity is proportional to the in¬ 
verse square of the radial distance. Our analytical and numerical solutions demonstrate 
decay of the orbits because of considering the drag force so that a cloud will eventually 
fall onto the central region after so called time-of-flight. We found that time-of-flight 
of a BLR cloud is proportional to the inverse of the dimensionless drag coefficient. We 
discuss if time-of-flight becomes shorter than the life time of the whole system, then 
existence of mechanisms for continually forming BLR clouds is needed. 
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1 INTRODUCTION 

In a unified theory for the structure of active galactic nu¬ 
clei (AGNs), various parts have been proposed to explain 
the observational features of t hese interesting astrophys- 
ical objects (e.g., iNetzerj I2OI3I) . Existence of cold clouds 
embedded in an intercloud medium known as broad-line 
region (BLR) may explain some of the observational evi¬ 
dences. Even near to the Galactic center, there are clouds 
in th eir orbits around the central supermassi ve black hole 
(e.g.. iGillessen et al.ll20l3 : iBurkert et al.ll2012l) . Apparently 
coexistence of clouds with a hot gaseous intercloud medium 
is common in some of the accreting systems. Current at¬ 
tempts to study these clumpy systems are concentrated on 
three main aspects. Understanding processes which may lead 
to th e formation of BLR clouds is an active resea rch field 
(e.g., lEromerth fc Meli^ 1200 ll : jPittard et al.l l2003l) . More¬ 
over, there are noticeable uncertainties about stability of 
these clouds and the confinement mechanisms in the light 
of theoretical arguments and recent numerical simula tions 
llReejflQ^ : iKrause et al]l2012l : fNamekata et al.ll2014l ). Re¬ 
gardless of current uncertainty over the nature of the con¬ 
finement mechanisms, orbital motion of the BLR clouds and 
their radiated emiss ion allow us to estimate mass of the cen¬ 
tral b lack hole (e.g.. [Marconi et al.ll200d : lNetzer &: Marzianil 
I2OIOI) . Considering different forces that may affect orbit of a 
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BLR cloud and the importance of the orbit’s shape, some au¬ 
thors addressed the o rbital analysis of the BLR clouds dur¬ 
ing recent years (e.g., [Krause et al.|[201ll : [Plewa et al.ll201^ : 
lKhaienabill2015l ). 

To make progress in this field, orbital motion of a cloud 
is treated like a classical two-body problem where gravita¬ 
tional force of the central black hole and a force due to the 
radiation of an accretion disc, as the dominant forces, are 
controlling t he net force on each cloud and the r e sultin g 
orbits (e.g., iLiu fc Zhangl [201II : [Krause et al.l [2OI1I . |2012|) . 
It is generally assumed that the clouds are optically thick. 
In another effort, just recently, a few authors studied dy¬ 
namical motion of clouds through a gaseous medium semi- 
analytically by assuming that the ensemble of clouds is col¬ 
lisionless and obtaining solution of the correspo nding Boltz¬ 
mann equation under sim plified conditions llWang et al.l 
I2OI2I : iKhaienabi et al.ll2014l) . Simplicity of treating the sys¬ 
tem as the two-body central force problem and its impli¬ 
cations on calculating the broad-line emission, however, de¬ 
serve further studies in analyzing orbital motion of BLR 
clouds in particular regarding to the unexplored physical 
ingredients. 

BLR clouds are subject to the anisotropic radiation 
pressure of a centr al source as has been pointed out by 
iLiu &: Zhangl ll201lj) . Recently, a few authors studied orbit 
of BLR clouds under combined effects of the central grav¬ 
ity and anisotrop i c radiation press ure (iKrause et al.l [201 ll : 
[Plewa et al.l l2013l : IKhaienabi I2OI5I ) . Eor pressure-confined 
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Figure 1. Location of a cloud is determined by the radial distance 
r and the polar angle 6. We assume the orbital plane of the cloud 
is defined by x-axis and OC so that the inclination angle is i and 

xdc = 'ip. 


BLR clouds with constant column density, IPlewa et alJ 
J2013t) presented analytical solutions for the orbit of the 
clouds. None of the previous studies of BLR cloud’s orbit 
considered interaction of a cloud with the ambient medium 
in the form of a resistive force. Although effect of the drag 
force on the orbital motion of the BLR clouds has not been 
studied, in the context of planetary motion and orbital anal¬ 
ysis of the satellites in a central gravitational force, a num¬ 
ber of authors treated the classical two-body problem with 


solutions fe.E., Humi & Carter! 20021: 

Mittleman & Jezewskil 

1 19821: iMavraeanis & MichalakisI 1994 

). 

iHumi & Carter ( 

2002|) studied 

properties of the or- 


bital motion of an object in a central force with a drag 
force opposite to the velocity vector and proportional to 
a power-law function of the magnitude of the velocity and 
the distance to the center of attraction. Based on their solu¬ 
tions, they demonstrated decay of the orbits when the drag 
force is considered. When the drag force is proportional to 
the velocity and the inverse square of the radial distance, 
iMittlem an fc Jezewskil l ll982fl obtain e d an analytical solu¬ 
tion. Then, Mavraganis fc Michalaki^ ll 1994 1 generalized the 
analysis to a case with an isotropic radiation force and found 
three integrals of motion. 

In this paper, we extend previous studies of BLR cloud’s 
orbital motion by considering a drag force which is opposite 
to the velocity of a cloud. We prescribe physical properties of 
the intercloud gas such as its pressure and dynamic viscosity 
as power-law functions of the radial distance. In the next 
section, general formulation of the problem is presented. In 
section 3, analytical solutions are presented for a given set 
of input parameters. Then, a detailed analysis of shape of 
orbits is performed in section 4. We conclude with possible 
astrophysical implications of our results in section 5. 


2 GENERAL FORMULATION 

We consider a pressure-confined cloud with mass m on its 
orbit in the gravitational field of a central object with mass 
M. The gas cloud is treated as a spherical solid body with a 
constant mass and no ablation. For simplicity, the intercloud 
gas is assumed to be static which means the rotational and 
the radial velocities of the gas flow are neglected. A cloud 


is subject to the other forces suc h as a non-isotropic radia¬ 
tion of the central accretion disc dLiu fc Zhangllioill l and a 
drag force in the opposite direction of its orbital motion. An¬ 
other reasonable assumption is that the volum e filling factor 
of the BLR cloud is low le.g.. lRees et al. which allows 

radiation pressure acts on each cloud directly and the shad¬ 
owing effects are neglected. Dependence of the drag force on 
the relative velocity of the cloud and the surrounding gas 
is assumed to be linear. This functional form of the drag 
force is known as Stokes’ drag which gives frictional force 
on a small spherical objects in a laminar fluid. Such kind 
of gas flow occurs at low velocities and the nondimensional 
Reynolds number is a diagnostic to determine if the flow 
tends to be laminar or turbulent. 

A primary uncertainty is that we do not yet know phys¬ 
ical properties of the intercloud gas of a BLR region. If we 
use Stokes’ drag relation, an estimate of the dynamics viscos¬ 
ity of the intercloud gas is needed. Although a BLR cloud is 
treated as a solid spherical object, it is actually a gaseous ob¬ 
ject which is possibly confined by an eternal agent like mag¬ 
netic field. Thus, it is not clear if dynamics viscosity which is 
used for the motion a solid body subject to the Stokes drag 
is adequate for the motion of a gaseous pressure-confined 
cloud through an intercloud medium. We can not resolve this 
uncertainty unless further numerical simulations for the mo¬ 
tion of a BLR cloud through an ambient gas are performed 
to obtain the drag force due to the interaction of the cloud 
with the intercloud. Thus, we approximate dynamic viscos¬ 
ity ^ using a standard relation, i.e. — (l/2)puA, where p, 
u and A are the density of the gas, the average molecular 
speed and the mean free path, respectively. Also, we have 
A ~ l/(d^p) where d is the typical diameter of an atom or 
molecule. On the other hand, Reynolds number is defined 
as Re = pvL/p, where v is the mean velocity of the object 
relative to the gas and L is a characteristic length. Thus, we 
have Re = 2{v/u){L/\). It is believed that the intercloud 
gas with the number density 10^cm“® is at a temperature 
aroun d or even larger than 10® K (e.g., IWhittle &: Saslaw 
Il986ll . So, we have u ~ 2 x 10®m/s and A « 10^°m. If we 
assume the mass of the central black hole is 10® Mq, then 
Keplerian velocity at 1 pc is approximately 6.5x 10"^m/s. Ob¬ 
viously, the velocity of a cloud may reach several 1000 km/s 
in the vicinity of the central black hole. But we approxi¬ 
mate the typical velocity of a cloud as u « 6.5 x lO^m/s. 
As for typical length scale, we approximate it as, say, ten 
times of the size of a BLR cloud, i.e. L ~ 10^®cm. Thus, 
we have Re ~ 0.65. Although our estimate of the Reynolds 
number provides a justification for the assumption that the 
drag force has a linear dependence on velocity, we have high 
velocity clouds at the inner parts of the system which imply 
a much larger Reynolds number and then a quadratic rela¬ 
tionship is more appropriate. But our analysis is still valid 
for the above mentioned conditions. Moreover, we did not 
consider the motion of the background gas. But if the BLR 
intercloud gas is rotating at a fraction of Keplerian velocity, 
then the relative velocity of a BLR cloud is reduced sig- 
nihcantly which implies a smaller Reynolds number. These 
issues should be studied in future studies. 

Figure [T] shows location of a cloud and our system of 
coordinates. Obviously, the angular momentum vector has 
a constant direction in space which implies the orbital plane 
has a fixed inclination angle i with respect to the accretion 
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disc, i.e. x — y plane. Thus, location of the cloud can also be 
written in terms of angle /XOC = ip in the plane of motion 
and the radial distance r. 

Thus, equations of the orbital motion become 

•2 GM GnyRci . 

r — rip = —— (/c sin'0 — 1)-r, (1) 

m 

rip + 2rip =(2) 

m 

where /i is the viscosity of intercloud gas and Rc\ is the radius 
of the cloud. Moreover, the non-dimensional parameter k is 
defined as 


k = 


31 

UmNcKJT 


sin(i). 


(3) 


Here, N^\ is the column density of a spherical cloud and /im is 
the mean molecular weight. The Eddington ratio is denoted 
by I = La/Tedd where La and Ledd are the luminosity of the 
central source and the Eddington luminosity, respectively. 
In the absence of drag force, the above equatio ns reduce to 
the main equations of the previous studies fe.g.. |Plewa et all 
I 2 OI 3 I : lKhaienabill2015ll . 

We note that in our model each cloud is in pressure equi¬ 
librium with the ambient gas. Then, radius of a cloud Rd 
becomes in proportion to Pgas where Pgas is the intercloud 
gas pressure. So, behavior of the pressure of the intercloud 
gas with radius is an essential part of the present model for 
the pressure-confined clouds. It is generally assumed that the 
intercloud pressure distribution is proportional to a power- 
law function of the radial distance as Pgas oc r~‘’ w here s 
is an input parameter fe.g.. lNetzer &: Marzianill2010l l. How¬ 
ever, this pressure profile may have a latitudinal depe ndence 
as it has been studied in detail bv iKhaienabil ll2015ll . Here, 
we do not consider such a dependence for simplicity. There¬ 
fore, radius of each cloud varies according to its location 
within the intercloud gas as 


/ ». \ s/3 

Rci = Rcio 

where ro is a reference radial distance, say, the initial radial 
distance of the cloud, and Lcio is the radius of the cloud 
at ro. For our pressure-confined clouds, the column density 
becomes oc or = No{r where No is a 
constant column density. As for the viscosity of the gaseous 
intercloud medium, we prescribe it as a power-law function 
of the radial distance, i.e. 


M = Mo 



(5) 


where u is an input parameter and yo is the viscos¬ 
ity of the intercloud gas at radius ro. Thus, our non- 
dimensional parameter k becomes k = ko{r/roY‘^^ where 
ko = 3lsin{i)/y^NoOT. It is assumed that the temperature 
of a clo ud does not chan ge significantly during its orbital 
motion llWang et al.ll201^ . 

Although our intention is not to study the confinement 
of BLR clouds by magnetic fields or other physical mech¬ 
anisms, there are some observational and theoretical moti¬ 
vations for the power-law functions which are assumed for 
the intercloud gas pressure and dynamic viscosity. Some au¬ 
thors also used similar power-law functions to investigate 
orbital motion of BLR clouds and obtained the emission 
profiles that are consistent with observed emission spec¬ 


tra le.g.. lNetzer fc Marzianill201Cll L Although our knowledge 
of physical properties of the intercloud gas is still poor, 
an idea put forward by a few authors that the intercloud 
gas can be described using Advection-Dominated Accretion 
Flows (ADAFs), in which physical qua ntities are power-law 
functions of the radial distance (e.g., iKrause et al.l I2OI2I : 
lKhaienabil2015ll . If the temperature of the intercloud gas 
is at the virial temperature, then the exponent of the pres¬ 
sure distribution is s = 5/2. But at the very inner parts of 
the system where the virial temperature is very high, the 
more appropriate value of the exponent is s = 3/2 because 
of the efficiency of cooling. Because of uncertainties in the 
values of the exponents, we will do a parameter study for 
different values of the exponents. We also found that for 
only one particular case with s = 0 and v = —2, there is an¬ 
alytical solution for the orbit of the BLR clouds. Although 
we do not have enough observational information to confirm 
these particular values, the analytical solutions give a bet¬ 
ter insight of what we may expect for other values of the 
exponents. 

In order to proceed, it is convenient to transform the 
main orbital equations into non-dimensional forms. In doing 
so, we use the reference radial distance ro to introduce a non- 
dimensional radial distance. Then, Keplerian velocity at this 
radial distance is written as rK(ro) = \JGM /ro and our 
time unit becomes to = ro/rK(ro). We now change variables 
as r —>■ ror and t —>■ tot. Thus, 

r — rip^ = — (^kor^"t^\ sin-i/)! — — ar'^'^^r, (6) 

and 

rip + 2rip =—ar'^'^'S'^^ip, ( 7 ) 


where a is a dimensionless drag coefficient, i.e. 

A 

a = 

Re 

Here, Re is Reynolds number of the initial condition and the 
non-dimensional parameter A is 


( 8 ) 


2^Rcio^>cio^^K(ro)^ 


(9) 


where uo is the typical velocity of the cloud and pcio is the 
initial density of the cloud. Moreover, the density of the 
atmosphere at the initial location of the cloud is denoted by 
Mo- 

Equations and © are our main equations for de¬ 
termining orbit of a cloud subject to a linear drag inside a 
gaseous medium with a power-law pressure distribution. R is 
unlikely to obtain analytical solutions for all possible values 
of s and iz. We can obtain illustrative analytical solutions, 
however, for certain values of these input parameters as we 
will show below. 


3 ANALYTICAL SOLUTIONS 

When variation of the intercloud gas pressure profile with 
the radial distance is not considered, we have s = 0. In the 
absence of the drag, t his particular case h as already been 
studied analytically bv IPlewa et al.l (|2013l L We now study 
this particular case, but considering the effects of the drag 
force. Thus, equation © can be written as 
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L + ar^'+^iP = 0, (10) 

where L = ip is the non-dimensional angular momentum 
per unit mass. If we set v = —2, this equation becomes 
integrable, i.e. 

L = Lq — aip, (11) 


where Lq is a constant. Thus, magnitude of the angular mo¬ 
mentum gradually decreases as a BLR cloud revolves around 
the central object. Total number of rotations before angular 
momentum becomes zero is Lo/(27tQ:). 

Under substitution u = 1/r and changing the indepen¬ 
dent variable from time t to angle ip, we transform our main 
equations © and 0 into a differential equation which is 
integrable. Obviously, we have 


dr 


1 du v! ■ 

= -2 V’: 

■'1^ 


( 12 ) 


dt V? dt 

where the prime denotes differentiation with respect to the 
angle ip. Thus, 


.. ;2 n” u' .. 

r^2 — ip - -Ip - -Ip. 

u-^ 

Orbital equations and 0 are re-written as 


(13) 


- —?/)- (-2 —Ip H- \-auip) -- 


= u{ko\sm‘ip\ - 1), 


(14) 


and 


— 2,^ip^ + — + awp = 0. (15) 

u 

According to equation (1151) . the coefficient of the ratio u'/u 
in equation dm vanishes. Thus, 

id 

m”-I-I t = —(1 — fcol sin'i^l). (16) 

ip^ 


We also showed that ip = u^(Lo 
equation becomes 

d^u 1 — feo| sini/)| 

dip'^ (Lo — aip)^ 


aip), and so the above 


(17) 


If we introduce a new variable 2 = (Lo — aip)ja, this equa¬ 
tion can be simplified further and it becomes 


1 - feo|sin(-^ - 2)1 


(18) 


If we neglect t he drag force, then e quation (1171) reduces 
to equation (7) in IPlewa et al.l ll2013f ) who studied orbits 
of BLR clouds without a resistive force. They obtained an 
analytical solution and a condition for having bound orbits. 
Their analytical solution is written as 


ko 

u(ip) = 1 -\- -^ip cos Ip -\- A sin ip + B cos i) 


(19) 


for 0 ^ Ip ^ n and similarly for tt < ip ^ 271 , we have 


u{ip) = 1 — —Ip cos ip + {A + ko) sin ip + {B + kon) cos ip, (20) 


where A and B are arbitrary constants and both u and 
du/dip ar e continuous a t ip — tt. Having the above analytical 
solutions, IPlewa et al. Ilioli) investigated orbits of the BLR 
clouds. Obviously, in the presence of drag force, the orbits 


would not be bound. Interestingly, we can find an analytical 
solution when the drag force is considered, too. 

Moreover, if we neglect non-isotropic force due to the 
centra l radiation field, equat ion jm reduces to equation 
(12) in iMittleman fc Jezewskil (Il982l l which is integrable and 
the solution is written in terms of the sine and cosine inte¬ 
grals. We can also present solutions of equation (1181) in terms 
of these integrals. For 2nn < ip < (2n -|- l)7r (n = 0,1, 2,...), 
the general solution of equation (O becomes 


m(2) = 

A 2 n cos( 2 ) -I- B 2 n sin( 2 ) — [sin( 2 )Si( 2 ) -I- cos( 2 )Ci( 2 )] 


ko 


cos(— -|- 2)Si(22) sin(— -|- 2)Ci(22) 
a a 


and for (2n -|- l)7r < ip < (2n -|- 2)7r, the solution is 


( 21 ) 


m(2) = 

C 2 n cos( 2 ) -I- D 2 n sin( 2 ) -|- ^ [sin( 2 )Si( 2 ) -I- cos( 2 )Ci( 2 )] 




cos(— -|- 2)Si(22) sin(— -|- 2)Ci(22) 
a a 


( 22 ) 


where A 2 n, B 2 n, C 2 n and D 2 n are arbitrary constants and 
the continuity of u and du/dz a± ip = 2n7i and ip = (2n -|- 
l)7r for a given n imply that not all these parameters are 
independent. 

Our analytical solutions clearly demonstrate decay of 
orbits because of the presence of the drag force. But shape 
of the orbits and the time for a BLR cloud to travel from an 
initial radius to the center depend on the initial conditions 
and the input parameters. 


4 ANALYSIS 

First, we study orbital shape of a BLR cloud in the plane of 
motion for a particular case with s = 0, u = —2, Lo = 1 and 
ko = 0.1 in Figure [2] We numerically solve the orbital equa¬ 
tions and our analytical solutions for this particular case 
are used to check the accuracy of our numerical solutions. 
Differences between our analytical and numerical solutions 
are found to be very negligible. In this figure and all sub¬ 
sequent figures, we consider the same initial conditions, i.e. 
r(t = 0) = 1, r(t = 0) = 0, ip{t = 0) = 0 and ip{t = 0) = 1. 
Moreover, the inclination angle of plane of motion i is fixed 
so that X axis coincides x axis and Y axis make angle i with 
y axis in y — z plane (Figure [T]). We also explored other 
initial conditions, but the results are qualitatively similar 
to what we are reporting here for a few illustrative cases. 
Figure [2] shows shape of the orbits for three values of the 
dimensionless resistive coefficient a. It is apparent that the 
decay of the orbit is much faster as a goes from 0.01 to 0.1. 
We can also calculate the time for a BLR cloud to travel 
from its initial location to the center which is defined as the 
time-of-fiight tf. Considering orbits in Figure [2] we can say 
that time-of-flight is significantly reduced with increasing 
the coefficient a. 

In Figure |3l the effect of the non-isotropic radiation 
force is studied for a configuration where s = 0, v = —2 
and a = 0.1. A case with ko — 0 corresponds to a system 
without a central source of the radiation. We see that the 
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X 



X 



X 

Figure 2. Orbit of a BLR cloud when its journey starts from 
the initial location r = 1 and the input parameters are s = 0, 
u = —2, kg = 0.1 and Lq = 1. Initial conditions are r{t = 0) = 0, 
i/;(i = 0) = 0 and = 0) = 1. Each plot is labeled by the 
corresponding value of a. 


orbits become more elongated as the parameter kg increases 
because of the non-isotropic nature of the central radiation 
field. Thus, time-of-flight tf increases as feo increases from 0 
to 0.6. 

We prescribed the background pressure as a simple 
power-law function of the radius with the exponent s. For 
our pressure-confined clouds, the exact mechanisms of the 



X 


Figure 3. Orbits of a BLR clouds when all the input parameters 
and the initial conditions are similar to Figure except that 
ko which varies from 0 to 0.6 in order to illustrate role of non¬ 
isotropic radiation in shape of the orbits. Each orbit is labeled by 
the corresponding value of ko. 



X 

Figure 4. Orbits of a BLR cloud for different values of s, i.e. 
s = 2.5 (solid-black), s = 1.5 (dashed-blue), s = 0 (dotted-red). 
The other input parameters are u = —2, a = 0.1, kg = 0.1, 
r{t = 0) = 0, ip{t = 0) = 0 and fit = 0) = 1. 


confinement are not considered as long as the pressure profile 
is prescribed as a power-law function of the radius. Figure 
|4] shows orbits of a BLR cloud for different values of the ex¬ 
ponent s while keeping all the rest of the input parameters 
unchanged. Solid, dashed and dotted lines are correspond¬ 
ing to s = 2.5,1.5 and 0, respectively. This Figure indicates 
that decay of the orbits depends on the exponent s so that 
time-of-flight tf becomes longer as the value of s increases 
from 0 to 2.5. Drag force is directly proportional to the ra¬ 
dius of a cloud which is a function of its location because 
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- 1.0 - 0.5 0.0 0.5 1.0 

X 


Figure 5. Orbit of a BLR cloud when the dimensionless drag 
coefficient is very small, i.e. a = 0.001. Input parameters are 
1 / = —2, /cq = 0.1 and s = 2.5 and the initial conditions are same 
as Figure |2] 


of pressure-confinement condition. Thus, drag force also de¬ 
pends on the radial distance so that a smaller value for s 
implies a larger drag force. For this reason, when we have 
s = 0, a BLR cloud very quickly falls onto the central object 
in comparison to the situations with the larger values for s. 

When the drag coefhcient is small, a cloud will decay 
very slowly but the non-isotropic radiation force causes the 
orientation of the orbit changes smoothly so that a cloud on 
its orbit may come back to same region for several times. 
In other words, presence’s probability of a cloud at specific 
regions is larger than other parts depending on the input 
parameters and under certain circumstances. As an illustra¬ 
tive example of this case, orbit of a BLR cloud for a small 
drag coefhcient but in the presence of non-isotropic radia¬ 
tion held is shown in Figure [5] We see that a cloud on its 
orbit repeats certain locations so that spiral patterns will 
emerge as a whole. If we start with an ensemble of identical 
clouds and the same initial conditions, the spatial distribu¬ 
tion of the clouds as they orbit will resemble to these spiral 
patterns. Although emergence of such a conhguration is ex¬ 
tremely unlikely and rather of academic interest, existence 
of these complicated orbits show richness and diversity of 
orbits when the drag force and the non-isotropic radiation 
held are taking into account in addition to the central grav¬ 
itational force. 

In our model, the viscosity is also prescribed as a power- 
law function of the radius with the exponent i/. We examine 
effect of varying parameter i/ on the shape of orbits in Fig- 
ure[6l Here, we assume that a = 0.1, ko = 0.1, and s = 2.5, 
but different values of u are considered. As the viscosity 
becomes steeper with the radius, then drag force becomes 
stronger which result in shorter time-of-hight. For other in¬ 
put parameters, dependence of the orbits on the variations 
of 1 / is found qualitatively similar to Figure [U] 

Prohle of t/ as a function a is shown for a particular 
case in Figure [T] The initial conditions are similar to the 


1.0 
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- 0.5 
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1 1 1 

- 1.0 

1 1 1 1 1 1 1 


- 1.0 - 0.5 0.0 0.5 1.0 

X 




X 


Figure 6. Orbit of a BLR cloud when we have ct = 0.1, fco = 0.1 
and s = 2.5. Each plot is labeled by the corresponding value of 
V. Again initial conditions are same as Figure [2] 


previous Figures and we have s = 0 and v = —2. Different 
values for feo are assumed. Curves from below to top are 
corresponding to feo = 0, 0.1, 0.2, and 0.3, respectively. This 
Figure shows that t/ oc a~^ and the time-of-hight becomes 
longer with increasing feo. 
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Figure 7. Time-of-flight tf as a function of the dimensionless 
drag coefficient a when the input parameters are s = 0 and u = 
—2 and the initial conditions are same as Figure [2] Curves from 
below to top are corresponding to kg equal to 0, 0.1, 0.2 and 0.3, 
respectively. 

5 DISCUSSION 

Our knowledge about origin of BLR clouds is not satisfac¬ 
tory. We do not know if these clouds are constantly created 
or they have formed more or less simultaneously some time 
in the past. In fact, there is an ongoing debate about if 
these clouds are stable structures or else they are transient 
objects which are constantly created and destroyed in a tur¬ 
bulent medium. It is well-known that magnetic fields play 
an important role in confining the clouds and if their effects 
are neglected the clouds would not be stable and very soon 
after their formation they will expand and disperse (e.g., 
Numerical simulations of clumps in a radiative 
medium under restrictive conditions such as neglecting mag¬ 
netic effects show that clouds are not able to move around 
because of the force du e to the radiation pres sure. On the 
other hand, just recentlv lMcCourt et al.l il2014l l showed that 
if magnetic fields and their role in confinement of the clouds 
are considered, the clouds not only become stable objects 
but also they can move in their orbits. As for the G2 cloud 
in the Galactic center, however, they predict that drag forces 
can lead to deviations from Keplerian profile. Thus, this is¬ 
sue deserves further investigations in particular regarding 
magnetic confinement and possible role of magnetic fields in 
modifying structure of the intercloud gas. But if the BLR 
clouds do really exist, they should be long-lived objects be¬ 
cause their dominant role in the emission from such systems. 
There is also a possibility that formation and destruction of 
the clouds are so fast that only on the basis of the current 
received emission we can not distinguish weather they are 
long-lived stable objects or not, unless more theoretical and 
observational diagnostics are proposed. 

Irrespective of the stability issue and the exact con¬ 
finement mechanism, we assume that the clouds are stable 
long-lived objects which have formed some time in the past. 
Because of the drag force, all clouds will eventually fall onto 
the central black hole or they may destroy at the very inner 


parts due to gravitational tide or intense central radiation. 
In any case, time-of-flight tf of a BLR cloud which starts 
its life at the outer parts of system should be longer than 
the life-time of the whole system, i.e. r. If the time-scale 
tf becomes shorter than the life-time of the whole system, 
one can expect all BLR cloud will eventually fall onto the 
central region in the absence of any continuous mechanism 
for creation of BLR clouds. In other words, the system will 
be depleted of BLR clouds in this case. But emission of sys¬ 
tem is interpreted in terms of existence of these cold clouds. 
So, there is a contradictory situation regarding to the exis¬ 
tence of BLR clouds if time-of-flight becomes shorter than 
the life-time of the system unless we propose robust mecha¬ 
nisms for continually creating BLR clouds so that even they 
fall onto the central object, the creation mechanisms lead 
to the formation of new clouds. The net effect is a clumpy 
system, though some of the clouds may be destroyed be¬ 
cause of falling onto central object. In our orbital analysis, 
time-of-flight tf is very sensitive to the input parameter, in 
particular dimensionless drag coefficient a has a vital role 
so that we can have very long or very short tf depending on 
this coefficient. 

If we set the life-time of the system, then tf should be 
longer than this life-time in order to have a clumpy system 
in the absence of creation mechanism. This constraint will 
gives us a constraint on the value of a which depends on local 
properties of intercloud gas. On the other hand, if we can 
constrain the input parameters including a independently, 
then we can compare tf to the life-time of the system. As we 
discussed above, a case with tf shorter than the life-time of 
the system, implies existence of mechanisms for continuous 
creation of BLR clouds. But if tf becomes longer than the 
life-time of the system such a formation mechanism may or 
may not exist. 

For a wide range of the input parameters, we found 
that tf ~ tfi(v,s)/a where the parameter to depends on 
the exponents v and s (Figure [T]). Now, if we have tf < r, 
then there might be a continuous mechanism for replenish¬ 
ment of BLR clouds. Thus, existence of such a mechanism 
implies a lower limit for the age of the whole system, i.e. 
Tmin — to{h',s)/a. If the age of system is older than Tmin, 
we can conclude that BLR clouds are continuously created 
and possibly destroyed. Using equations (|H]| and ([21 , we then 
have Tmin = {2Re/9)toR, s)(Rcio/ro)(pcio/po) where we have 
assumed that uq ~ UK(ro). Observations suggest a BLR re¬ 
gion has a typical size of ~ 0.01 — Ipc which contains a large 
number of clouds with a very low filling factor (~ 10“®) 
and typical sizes ~ 10^^ — 10^"*cm and mags ~ 10~®M( 7) 
dRees et al.l Il989l : iMarconi et al.l l2008l : iPIewa et al.l l2013ll . 
Now, we can estimate Tmin if we set Re ~ 1 and assume that 
a cloud starts its journey from the outer parts of the BLR 
region, i.e. tq = O.Olpc. Also, we have Rdo = 10^"*cm. Thus, 
we obtain Tmin = 7.2x 10“'*(pcio/po)lo(i^, s). A larger density 
contrast Pdo/po implies a smaller cloud if the mass is kept 
fixed. If we set pdo — 10^°cm~® and po — 10‘^cm“® (e.g., 
[whittle fc Sasla^ ll986l L then we have Tmin = 720to(*^, s). 
Although value of to(iz, s) depends on the exponents u and 
s, its variation with these parameters is not very significant 
so that we generally found that toR, s) < 1. Thus, if we as¬ 
sume to{v,s) ~ 1, then Tmin = 720 or Tmin — 15 X 10®yr for 
a cloudy BLR region around a central black hole with mass 
10 ®Mq. But it is unlikely that a BLR region to be younger 
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than this time scale. In estimating Tmin, we considered the 
input parameters so that the largest possible value for this 
time-scale is obtained. In other words, the actual value of 
Tmin is probably much shorter than 15 million years unless 
the clouds become very compact with a very large density 
contrast which is not supported by observations. 

Line profiles of the BLR clouds based on their orbital 
motion has been calcul ated by some authors in the absence 
of the drag force le.g.. l^tzer fc Marziani| 120101 1. Consider¬ 
ing the importance of the drag force in changing the orbits, 
as we showed in this study, future work is also needed to 
calculate line profiles of the BLR clouds according to our 
orbital solutions. Moreover, we assumed that the drag force 
is a linear function of the velocity. Although it enabled us 
to obtain analytical solution for a certain set of the input 
parameters, the next step is to obtain orbits of BLR clouds 
when drag force is proportional to the velocity square of a 
cloud. 


6 CONCLUSION 

We investigated orbital motion of the BLR clouds which are 
subject to the gravitational force of the central black hole, 
a non-isotropic force due to an accretion disc and a resis¬ 
tive force in proportion to the relative velocity of the cloud 
and the intercloud gas. Although direction of the angular 
momentum vector does not change, its magnitude reduces 
because of the effect of the drag force and therefore this force 
inevitably causes the clouds to spiral downward. Assuming 
that the physical properties of the intercloud medium like 
its pressure and the viscosity are described as power-law 
functions of the radial distance, we found that the orbit of 
a cloud and its time-of-flight sensitively depend on the non- 
dimensional parameter a which is a function of the prop¬ 
erties of gaseous component at the initial location of the 
cloud. Our main conclusion is that most of the BLR clouds 
will fall onto the central black hole because their time-of- 
flight is shorter than life-time of whole system. This result 
is more or less independent of the exponents of the pressure 
distribution and dynamic viscosity. A constant replenish¬ 
ment mechanism for the BLR clouds is needed when the 
drag force is considered. 


ACKNOWLEDGMENTS 

This work has been supported financially by Research Insti¬ 
tute for Astronomy & Astrophysics of Maragha (RIAAM) 
under research project No. 1/3720-60. I am very grateful to 
referee for his/her report which greatly improved quality of 
the paper. 


REFERENCES 

Burkert A., Schartmann M., Alig C., Gillessen S., Genzel 
R., Fritz T. K., Eisenhauer F., 2012, ApJ, 750, 58 
Fromerth M. J., Melia F., 2001, ApJ, 549, 205 
Gillessen S., Genzel R., Fritz T. K., Quataert E., Alig G., 
Burkert A., Guadra J., Eisenhauer F., Pfuhl O., Dodds- 
Eden K., Gammie G. F., Ott T., 2012, Nature, 481, 51 


Humi M., Garter T., 2002, Celestial Mechanics and Dy¬ 
namical Astronomy, 84, 245 
Khajenabi F., 2015, MNRAS, 446, 1848 
Khajenabi F., Rahmani M., Abbassi S., 2014, MNRAS, 
439, 2468 

Krause M., Burkert A., Schartmann M., 2011, MNRAS, 
411, 550 

Krause M., Schartmann M., Burkert A., 2012, MNRAS, 
425, 3172 

Liu Y., Zhang S. N., 2011, ApJL, 728, L44 
Marconi A., Axon D. J., Maiolino R., Nagao T., Pastorini 
G., Pietrini P., Robinson A., Torricelli G., 2008, ApJ, 678, 
693 

Mavraganis A. G., Michalakis D. G., 1994, Celestial Me¬ 
chanics and Dynamical Astronomy, 58, 393 
McCourt M., O’Leary R. M., Madigan A.-M., Quataert E., 
2014, ArXiv e-prints 

Mittleman D., Jezewski D., 1982, Celestial Mechanics, 28, 
401 

Namekata D., Umemura M., Hasegawa K., 2014, MNRAS, 
443, 2018 

Netzer H., 2013, The Physics and Evolution of Active 
Galactic Nuclei 

Netzer H., Marziani P., 2010, ApJ, 724, 318 
Pittard J. M., Dyson J. E., Falle S. A. E. G., Hartquist 
T. W., 2003, A& A, 408, 79 

Plewa P. M., Schartmann M., Burkert A., 2013, MNRAS, 
431, L127 

Rees M. J., 1987, MNRAS, 228, 47P 
Rees M. J., Netzer H., Ferland G. J., 1989, ApJ, 347, 640 
Wang J.-M., Cheng C., Li Y.-R., 2012, ApJ, 748, 147 
Whittle M., Saslaw W. C., 1986, ApJ, 310, 104 

This paper has been typeset from a Tj/X/ MjjiX file prepared 
by the author. 


© 2015 RAS, MNRAS 000, [HS] 




